Non-equilibrium behavior of lysozyme solutions: beads, clusters and gels 
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Observation of salt-induced aggregation of lysozyme at pH = 4.5, 22° C by optical microscopy 
revealed four regimes: bicontinuous texture, 'beads', large aggregates, and transient gelation. The 
interaction of a metastable liquid-liquid binodal and an ergodic to non-ergodic transition boundary 
inside the equilibrium crystallization region can explain our findings. Lysozyme data at -2°C, as 
well as many literature results, are consistent with the relative movement of these boundaries. 



Short-range interparticle attractions are ubiquitous in 
synthetic colloids. Globular proteins in salt solution can 
also be modelled as 'sticky particles' 1, 2, 3]. The equilib- 
rium phase behavior of particles with an isotropic short- 
range attraction is well known: the attractive interaction 
widens the fluid-crystal coexistence gap. When the inter- 
particle potential is quantified by the second virial coef- 
ficient, a quasi-universal fluid-crystal coexistence bound- 
ary exists for synthetic colloids and proteins 4 S B| • 

The non-equilibrium behavior of sticky particles is far 
less understood. In the case of colloids, it is known that 
two qualitatively different glassy states exist at volume 
fractions ^ > 0.6 7]. When ^ < 0.2 and the attrac- 
tion is deep enough, permanent gelation is almost in- 
variably observed. At shallower attractions (< IO/cbT), 
'transient' gelation and separation into particle-rich and 
particle-sparse amorphous phases have also been de- 
scribed 8, 9, 10], together with a recent report of a 'clus- 
ter phase' Similar behavior occurs in protein solu- 
tions, although gels are less commonly reported. The- 
oretically, the gelation of sticky particles has long been 
discussed in terms of diffusion-limited cluster aggregation 
(DLCA) and percolation; recently mode-coupling theory 
(MOT) has been brought to bear on the problem . It 
is also known that a metastable fluid-fluid phase bound- 
ary (binodal) exists. In polymer solutions, the interfer- 
ence of gelation and fluid-fluid phase separation leads to 
a rich 'zoo' of non-equilibrium behavior ^^^^ 
Letter, we show that this is also the case with particles. 

Our particles are the globular protein, lysozyme. 
Sodium chloride induces an effective short-range attrac- 
tion by screening the intermolecular electrostatic repul- 
sion Careful visual observations and microscopy 
delineate four non-equilibrium regimes. We give a unified 
description of these observations within a framework that 
should be applicable to other sticky-particle systems. We 
show that the precise kind of phenomena observable in a 
particular system is sensitive to a range of physical pa- 
rameters. This helps explain why a plethora of reported 
observations exists on apparently similar systems. 

Six-times crystallized chicken egg-white lysozyme 
(Seikagaku America) was dissolved directly in a 50mM 
sodium acetate buffer titrated with hydrochloric acid to 
pH = 4.5. Stock solutions at 150 mg/ml were centrifuged 
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FIG. 1: Experimental phase diagram of lysozyme and NaCl at 
22° C and pH = 4.5. Axes give salt concentration, Cs, in mo- 
lars, and lysozyme concentration in volume fraction, cj) and 
weight /volume, c. Homogeneous solution (o); crystals (x); 
bicontinuous texture (■); spherical beads (a); large aggre- 
gates (a); transient gels (*). Continuous lines estimate the 
positions of the crystallization boundary lowest), non- 

equilibrium boundary (NEB, middle) and gelation boundary 
(upper, piecewise straight), the low-cj) portion of which ex- 
trapolates to ~ 0.13 (thin line). Inset: The height h (scaled 
to 1 at t = 0) of a transient gel (0 = 0.07, = 1.4 M) vs. 
time, t (in minutes). The 'latency time' is ~ 30 mins. 



to remove dust and undissolved material. Higher concen- 
trations were reached using a Vivaspin 6 ml concentrator. 
Exact concentrations were determined by UV absorption 
spectroscopy after 1000 x dilution in buffer, using a spe- 
cific absorption coefficient of 2.64 ml/mg.cm. (We cal- 
culate an effective (/) using a molecular radius of 1.7 nm 
and molecular weight of 14320 g/mol.) Salt solution fol- 
lowed by deionized water were pipetted to the protein in 
buffer, mixed by shaking and observed by the naked eye. 
Small amounts were drawn into O.lmm-thick rectangular 
capillaries and observed using a Zeiss optical microscope. 

At increasing salt concentration, c^, samples remained 
homogeneous (o, Fig.Q) until crystals (x, Fig.^ nucle- 
ated beyond a crystallization boundary (CB). At higher 



Cg, across a 'non-equilibrium boundary' (NEB), all sam- 
ples turned turbid after mixing. In the microscope, four 
regimes can be distinguished (■, A, A and *, Fig.Q). 

Macroscopically, samples with the highest Cg (*, Fig.0 
turned white throughout after mixing. A plot of the 
height of the white portion of the sample vs. time showed 
a reverse sigmoidal shape (inset, Fig.P). In other words, 
rapid sedimentation would start suddenly after a 'latency 
period' (~ 15-30 minutes). Fig. inset. The result- 
ing sediment crystallized after a few hours. Such 'de- 
layed sedimentation' is characteristic of 'transient gels' 
P . Note that our lysozyme transient gels are inhomoge- 
neous on the ~ 10/im (~ 10^-particle) level, Fig.[2fe). 

Non-gelling samples above the NEB turned cloudy af- 
ter mixing, and gave rise to turbid sediments without 
delay times. Although to the naked eye there is little to 
distinguish between the different samples in this region, 
microscopy revealed three distinct regimes. Just across 
the NEB and at (/) > 0.1 (■, Fig.[H), we observed bicon- 
tinuous textures reminiscent of spinodal decomposition. 
Fig. El^a). Crystals subsequently nucleated exclusively 
in one of bicontinuous regions. Samples just across the 
NEB but with (j) < 0.1 (A, Fig.Q) showed freely- moving, 
non-coalescing spherical 'beads' a few jam in diameter, 
or small aggregates of these beads. Fig. |2jb). Finally, at 
higher Cg (A, Fig.CJ, we observed progressively larger ag- 
gregates; within each aggregate, structures on the scale 
of the previously-mentioned 'beads' could be discerned, 
Fig|2fc). The beads and large aggregates would slowly 
dissolve once crystallization started within a few hours. 

The interaction between lysozyme molecules at high Cg 
can be modelled by an isotropic, short-range attraction 
Q, 1^, m . Our observation of gelation at the highest Cg 
(*, Fig.^, which can be compared to the lysozyme gels 
observed by Muschol and Rosenberger 0, is therefore 
unsurprising. In an idealized system of particles with 
infinitely short-range and infinitely deep interparticle at- 
traction, we expect irreversible DLCA to lead to a space- 
spanning gel jla,ll3 • Gelation in real systems with finite- 
range and finite-depth attractions is well known, but its 
cause is intensely debated 18j, a possible explanation 
being an ergodic to non-ergodic transition predicted by 
MCT ^3 or a 'renormalized' variant jl9| . 

The non-permanent interparticle bonds give rise to 
ageing and, ultimately, collapse of the gel structure and 
rapid sedimentation after a 'latency period'. Such 'de- 
layed sedimentation' is ubiquitous in synthetic colloids 
with short-range attraction (see j3, U and references 
therein). After our lysozyme gels collapsed, crystal- 
lization occurred. 'Transient gelation' is therefore a 
metastable state en route to thermodynamic equilibrium. 

The bicontinuous spinodal- type textures (■, Fig.^and 
Et) are also unsurprising. There is a metastable fiuid- 
fiuid binodal 'buried' within the equilibrium fluid-crystal 
coexistence boundary of sticky particle systems 0, , 
including lysozyme in salt solutions 
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FIG. 2: Phase contrast images (symbols as in Fig.Q) showing 
(a) bicontinuous texture, (0 = 0.14, Cs = 0.7 M); (b) spheri- 
cal beads (note higher magnification), (0 = 0.07, Cs = 0.9 M); 
(c) large aggregates, (0 = 0.02, Cs = 1.5 M); (d) large aggre- 
gates and a crystal (bright, middle), (0 = 0.05, Cs = 1.1 M), 
between crossed polarizers; (e) transient gel, (0 = 0.1, Cs = 
0.9 M), with inhomogeneities on the ^ lO/xm scale. 



The 'beads' and large aggregates are more intriguing. 
NMR measurements of these detected two populations 
of proteins: a fraction of free molecules (presumably in 
solution), and a fraction with severely restricted molec- 
ular motion. Preliminary X-ray diffraction gave no pow- 
der rings. Under crossed polarizers, the beads and large 
aggregates appeared birefringent, although less so than 
crystals, Fig|2{d). These results suggest that the beads 
and the large aggregates consist of lysozyme in an amor- 
phous solid state. Our 'beads' are reminiscent of the 
'oil-like droplets' frequently reported in studies of protein 
crystallization. A recent review 22] noted that 'crystals 
... nucleate on the surfaces of [these] droplets . . . and they 
always grow outward into the exterior . . . protein contain- 
ing medium, not inward into the . . .protein rich phase.' 
This observation suggests that these 'oil droplets' may in 
fact be amorphous solid phases on which heterogeneous 
nucleation takes place, consistent with our 'beads'. 

We suggest that the formation of 'beads' and large ag- 
gregates is due to the interaction of kinetic arrest and 
liquid-liquid phase separation. Our proposal is summa- 
rized in Fig.lSl Here CB is the equilibrium crystallization 
boundary and LL is a liquid-liquid binodal. An ergodic 
to non-ergodic transition ('kinetic arrest') occurs across 
NErg, i.e., this is where gelation would occur in the ab- 
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FIG. 3: Schematic illustration (in color on line) of the ori- 
gins of observed nonequilibrium behavior for samples ob- 
tained by diluting protein solution with volume fraction A 
with salt solution. CB — crystallization boundary, NErg = 
(underlying) ergodic to non-ergodic transition boundary, LL 
= (metastable) liquid-liquid binodal with critical point (CP) 
and tie lines ab, cd. We predict gelation above dg-gC\ large 
aggregates and beads in the region between gC' and gc, liquid- 
liquid phase separation below dc but above LL, and crystal- 
lizaton between LL and CB. 



sence of LL. The nature of this transition is debated 
Ellillii- At ^ > 0.4, its position may be given rather 
accurately by MCT 0. However, the precise identity of 
NErg does not affect what follows. 

Now recall that a sample with (average) composition 
e would have been obtained by mixing protein solution 
A and brine A'. We suggest that in this process the 
protein solution is 'quenched' along the path AA' 0. 
When the quench path crosses LL at a, phase separation 
into protein-rich and protein-poor solutions with compo- 
sitions given by the ends of the tie line ab occurs. If the 
quench is sufficiently slow compared to the rate of this 
phase separation, there will be time for macroscopic do- 
mains to develop and densify following the binodal along 
ad and be. In other words, a slow quench permits us to 
assume that a sample is always in local thermodynamic 
equilibrium. If the average composition is somewhere on 
ae, continuously coarsening bicontinuous texture results, 
until one of the phases crystallizes (cf. ■, Fig.Q]). 

If the average composition is e, the bicontinuous 
texture coarsens until the compositions of the phase- 
separated regions reach c and d. Here, the continu- 
ous protein-rich region will become non-ergodic: d is on 
NErg^ and we find a gel with spatial inhomogeneities gov- 
erned by the length scale of the bicontinuous texture at 
the point of gelation (cf. *, Fig. ^ and Fig.I^fe)). 

Along the tie line dc there exists a 'dynamic perco- 
lation threshold' 24], point to the left of which the 
phase-separated texture is no longer bicontinuous. Thus, 



for a quench path just to the left of AgC\ we expect 
a texture that just fails to be connected through space. 
When the protein-rich portions of this locally bicontin- 
uous texture becomes non-ergodic, large aggregates re- 
sult (A, Fig. d. Paths further to the left of AgC give 
more disconnected phase-separated domains and there- 
fore smaller non-ergodic aggregates. Thus, quenching 
along a path from A towards / should give more or less 
isolated domains: protein-rich droplets. Their Ostwald 
ripening gives a peaked size distribution 25]. When the 
composition of these droplets reaches d, they turn non- 
ergodic to give rather monodispersed 'beads' (A, Fig.QJ, 
whose mean size is determined by the interplay between 
the quench speed and the phase transition kinetics. 

Thus, our prediction is that for all samples prepared 
by diluting protein solution A with salt solution A', and 
at protein concentration below d^ dg-gC is the gelation 
boundary with gC extrapolating back to A on the pro- 
tein axis, and LL is the NEB. Experimentally, our gel 
boundary (separating * from other symbols in Fig. is 
consistent with this prediction: it is plausibly piecewise 
straight. The observed NEB (the upper boundary of x. 
Fig. ^ indeed has the shape of a LL binodal. 

In detail, all our samples with ^ < 0.1 were prepared 
by diluting protein solution of volume fraction = 0.132 
with brine. Consistent with Fig. [31 the steeper part of the 
experimental gel boundary extrapolates to (j) = 0.13 ~ 
00- Again, according to Fig. [3 the shallower part of the 
gel boundary is part of a tie line of the LL binodal. If 
we assume this interpretation for the experimental gel 
boundary in Fig. Q then two deductions follow. First, 
we can estimate (using the lever rule) that the sample at 
(/) = 0.14, C5 = 0.7 M (right-hand ■) should give - 2/3 
high-protein-concentration domains, consistent with mi- 
croscopic observations. Fig. [2fa). Secondly, we can esti- 
mate that the critical density 0crit ~ 0.12 — 0.13, consis- 
tent with (/)crit ~ 0.15 under similar conditions (151] . 

Finally, quench-path dependence is a crucial part of 
our explanatory scheme: the observed gel boundary 
should change if the starting protein concentration (and 
therefore the quench path) is altered. A sample that 
gave large aggregates, [a]. Fig. ^ gelled when prepared 
by diluting with brine a protein solution at = 0.176. 

To summarize, we have observed different regimes of 
salt-induced aggregation and gelation in lysozyme so- 
lutions at pH = 4.5 and 22° C. The existence of these 
regimes can be explained by a metastable liquid-liquid 
binodal and a non-ergodicity transition line intersecting 
within the equilibrium fluid-crystal coexistence region. 

It is probable that the relative positions of the vari- 
ous boundaries shown in Fig. 01 and therefore the kinds 
of non-equilibrium behavior observable in an experiment, 
are rather sensitive to the precise parameters of a particu- 
lar system (e.g. the range of the attraction 113), which in 
turn should be temperature dependent 01 • We repeated 
our experiments at temperatures -2°C to 40° C. Details 
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FIG. 4: Experimental phase diagram of lysozyme and NaCl 
at — 2°C and pH = 4.5: crystals (x) and gels (*). 



will be reported elsewhere . Here we show results at - 
2°C, Fig.^ where the NErg boundary has dropped below 
the LL binodal. Gelation was the only non-equilibrium 
behavior observed. 

We believe that the scenario described in Fig.|31is also 
relevant to other experiments 0, [2^ and simula- 
tions reported 2"^ in the literature . In many exper- 
iments, however, there are not enough published data 
points to justify this claim to the level of detail possible 
using the data in Fig. ^ More importantly, many data 
sets pertained to colloids with diameter in the 0.1 — 1/im 
range. The larger size slows down the rate of phase sep- 
aration relative to the intrinsic particle dynamics (thus 
giving slower quenches). It also increases strongly the 
effect of gravity: clusters may sediment before they have 
time to grow large enough and touch to form a system- 
spanning gel. In some cases, charge may also be present, 
and the balance between a very short-range attraction 
and a longer range repulsion could generate non-trivial, 
new phenomena 30]. A detailed study of charge and 
gravity effects will be presented elsewhere 31j. 

Thus, while there may be universality in the equilib- 
rium phase behavior of particles with short-range attrac- 
tion 0, H, 0, 5 we have demonstrated that a bewil- 
dering 'zoo' of non-equilibrium phenomena is possible, 
depending on precise conditions. A scheme involving ki- 
netic pathways governed by the interaction between ther- 
modynamic and kinetic transitions has been proposed 
that is able to capture and conceptualize this complex- 
ity. 

We thank Lindsay Sawyer for X-ray diffraction and 
Mike Gates for discussions. KK was funded by the EU. 
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